The quantum valley Hall effect in proximity-induced superconducting graphene: an 
experimental window for deconfined quantum criticality 
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We argue that by inducing superconductivity in graphene via the proximity effect, it is possible 
to observes the "quantum valley Hall effect". In the presence of magnetic field, supercurrent causes 
"valley pseudospin" to accumulate at the edges of the superconducting strip in Fig. This, and the 
structure of the superconducting vortex core, provide possibilities to experimentally observe aspects 
of the "deconfined quantum criticality" . 

PACS numbers: 



Neutral graphene is a semimetal with Fermi points 
(K + and K_) at the Brillouin zone corners ||. Van- 
ishing density of states and low dimensionality tend to 
prevent graphene from developing order. However, it 
has been demonstrated that through proximity effect by 
making a contact with SC electrodes, graphene can carry 
dissipationless supercurrent || Q which is a promising 
observation for possibility of inducing superconductivity 
through proximity effect in graphene|5| (Fig. [l]). 

In this letter we show that if one views the valley in- 
dex K± as a pseudospin degree of freedom (valley pseu- 
dospin), the proximity induced SC graphene will exhibit 
the "quantum valley Hall effect" (QVHE) in the presence 
of magnetic field. The QVHE is a close cousin of the 
more familiar quantum spin Hall effect (QSHE) J(| 
Here, supercurrent induces valley spin accumulation at 
the edges of superconducting graphene. 

This effect is closely related to the physics of the so- 
called deconfined quantum criticality^) where continuous 
Landau-forbidden transition is made possible because the 
defects of one phase carry the quantum number of the 
order parameter of the other phase. Interestingly de- 
spite the fact that neither superconductivity nor "Kekule 
distortion" is realized (or even nearly being realized) 
in graphene, the vortex core in the proximity induced SC 
graphene carry Kekule texture under suitable condition 
(see later). This "duality" between the superconduct- 
ing and Kekule order is in sharp contrast with the usual 
competing order phenomena where nearly degenerate or- 
ders compete for realization in the ground state. As a 
consequence when one order is suppressed in the vortex 
core the other order emerges [fL0| . If the SC order and the 
Kekule distortion are somehow stabilized in graphene, 
the above duality will make the continuous transition be- 
tween the Kekule order and superconductivity possible. 

The physical system we consider consists of a bulk s- 
wave superconductor coated by a single layer graphene 
shown in Fig. |l|. Through the proximity effect, SC can be 
induced in graphene Q. In addition since quasiparticles 
in graphene carry multiple quantum numbers, different 
types of paring could be induced. In the following we 




FIG. f : Proximity induced superconductivity in graphene. 
(a) Quantum valley Hall effect and (b) superconducting vor- 
tex. In (b), small staggered potential is introduced near the 
core to detect accumulated valley spin. 



consider the simplest on-site s-wave pairing described by 
the following Hamiltonian 



H = 
h 

A*A~ 1 



AA 

-h T 



(l) 



Here ip s = (ipA-ts,ipB+s,i>A-s,'>pB-s) is the four com- 
ponent fermionic field with spin s. ipA/B,+/-.s represents 
the electron annihilation operator on sublattice A/B as- 
sociated with valley K± and spin s. In the following, 
we use three sets of Pauli matrices, {r x ,y,zfl}, Wx^.z.o}, 
and {p x ,y,z,o}, each acting on sublattice (A,B), valley 
(+, — ) and Nambu (ip^, V'{ T ) degrees of freedom, respec- 
tively. In Eq. (|l|), ft, is given by h = diag (h+, h-) where 
h± = —i(r x d x ±Tydy), with the Fermi velocity set to one 
for convenience; the term A A with A G C and A = cr x r z , 
pairs electrons associated with K + and K_. 

Now wc examine the core structure of a superconduct- 
ing vortex in the SC graphene. The vortex is created by 
applying magnetic field perpendicular to the graphene 
plane, which induces vortex lines in the bulk SC and 
pass through the graphene layer. With a single vortex 
(with vorticity n) centered at r = the A in Eq. (pi) 
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becomes A(|r|)e where A(|r|) oc |r|" as |r| — > and 
A(|r|) — > A as |r| — > oo. In addition the magnetic field 
is introduced via the vector potential through minimal 
coupling. In the presence of the vortex the BdG equa- 
tion K^>(r) = E^!(r) acquires two zero energy solutions 
(for simplicity we consider single winding vortex) |ll| . In 
the symmetric gauge A(r) = -A(|r|) eg, the zero energy 
eigenfunctions ^lfr) and ^ 2 (r) are given by: 

ff(r) = [/(|r|),0,0 ) 0,0,0 J 0,-i/(|r|)], 

= [0,0,0, -i/(|r|),/(|r|), 0,0,0], (2) 

where the localized function /(|r|) is given by /(|r|) = 
e iSo A ( s ) ds g(\r\) with g(|r|) — ► constant as |r| — > and 
g(\r\) — > e~ A ' r ' as |r| — > oo (q > is the electric charge). 

Interestingly these degenerate vortex core states form 
the invariant space of a valley spin SU(2) symmetry as- 
sociated with the valley degeneracy of the BdG Hamilto- 
nian. The three generators of this SU(2) are 

Qi=T x a y p , Q 2 = -r x a x p z , Qz = T a a z p z , (3) 

with the expectation value of Q3 equal to the difference 
of the electron number associated with the two different 
valleys K±. It is easy to check that [Qi,Qj] = 2ieijkQj- 
Qz is diagonal in the representation given by Eq. (||) , i.e. 
the two zero modes carry ±1 Q3 quantum numbers. In 
the presence of valley spin SU(2) breaking perturbation 
the zero-mode degeneracy will be lifted. Within the sub- 
space spanned by the zero modes such effect is described 
by an effective magnetic field, B e s ■ T, where three com- 
ponents of r are the valley spin Pauli matrices. With 
B c g ^ the SU(2) symmetry is broken down to U(l) 
and we expect the split zero modes to be eigenstates of 
the generator of the residual U(l). 

The above SC vortex core physics is nicely captured 
by the following double Chern-Simons action 

S C s[A, A 5 , 6] = ^ / dtdre^Ald, (d p Q + 2qA p ) , (4) 

where A is the ordinary (external) electromagnetic gauge 
field, O is the phase of the SC order parameter, and 
A$ is the gauge field that minimally couples to the 
residual U(l) phase of the valley spin. By taking the 
derivative with respect to the temporal component of 
A\ in Eq. (§), we obtain AQ 5 = 5S C s/SA%\ A 5 =0 = 
(2tt)^ 1 J dre i: > didjO = n, where n is the vorticity, and Q 5 
is the generator of the residual U(l) valley spin symme- 
try. Formally, the action (||) can be obtained by integrat- 
ing out fermions in the presence of arbitrary background 
of A, A 5 and 6 (see Ref. @). The double Chern-Simon 
term (Q) is the formal declaration of the quantum valley 
Hall effect discussed in the introduction. 

With valley spin SU(2) symmetry breaking the SC vor- 
tex core acquires an interesting texture which allows its 
SU(2) structure to be identified by experiments. As we 
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TABLE I: Valley pseudospin and electric charge generators 
(left column) with corresponding order parameters (Higgs 
mass terms) which are rotated by the generator (right col- 
umn). 

shall discuss below, this texture is in the form of site or 
bond charge density wave (CDW). The site CDW breaks 
the equivalence between the A and B sublatticc without 
enlarging the unit cell. It is described by a real order 
parameter. On the other hand the bond density wave 
(or Kekule order|J) triples the unit cell size, and is de- 
scribed by two real order parameters each correspond- 
ing to one of the two complementary Kekule distortion 
patterns^. In the presence of these three types of order, 
the following "mass" terms / rfr$*M„$ are added to the 
Dirac Hamiltonian, where Mcdw = M3 = o~QT z p z , and 
M Ke ki = Mi = aiT 2 p 0l M Kck2 = M 2 = o 2 r 2 p z . 

It turns out that each of the Qi.,2,3 generates a rotation 
of a pair of the density wave order parameters discussed 
above, just like the electric charge operator Q generates 
a rotation between the real and imaginary part Mr c sc = 
o~ x T z p x , Mi m sc = o~xT z p y of the SC order parameter (see 
Table g) . 

One can infer these relations by noting that under a 
uniform rotation around i-th axis Qi in the valley spin 
space VP — > e* ¥ ' i< * i ^ r , the Hamiltonian in the presence of 
Higgs mass terms Mj k [where k) is a cyclic permu- 
tation of (1, 2,3)] (i.e. H + J dr^ (rijMj + n k M k ) *) is 
invariant if the rotation in ^ is accompanied by a simul- 
taneous rotation in rij t k, Tij + in k — > e~~ 2%Vi (nj + ink) i 
as one can check from the commutation relations, 
[Qi, Mj] = 2ieijkMk. In turn, this implies the presence of 
valley spin SU(2) breaking perturbation induces a linear 
combination of the three types of order discussed above. 
Within the subspace spanned by the vortex zero modes, 
it turns out that B e ff,i,2,3 = "1,2,3- 

We have studied two different types of valley SU(2) 
breaking by solving the BdG Hamiltonian numerically on 
a finite graphene lattice. In the first case, the interaction 
of the vortex with the edge of the system induces a B e s 
pointing in the x-y plane leading to a Kekule type vortex 
texture. In Fig. |3j we plot the expectation value of the 
nearest neighbor link operator cj s c- s in the core states 

(cj s is the lattice electron operator at i-th site with spin 
s). Here blue (red) denote positive (negative) values and 
the thickness represent the amplitude. It is interesting to 
note that the Kekule texture associated with the two split 
core states are opposite. The second type of valley SU(2) 




(a) 



(b) 



FIG. 2: Kekule texture observed for the split core states at 
(a) positive and (b) negative energy. Blue (red) bonds show 
strengthened (weakened) bonds, respectively. The filled circle 
denotes the position of the vortex. 



breaking is triggered by explicitly turning on a staggered 
chemical potential. This can be achieved through the in- 
teraction with the "buffer layer" when graphcne is grown 
on a substrate(l3| (Fig. ^(b)). In Figs, ga] and |3(b)| we 
show the induced site CDW texture associated with the 
two core states. Similar to the Kekule texture discussed 
above, the site CDW texture associated with the two core 
states are opposite. In principle these vortex core texture 
can be probed by scanning tunneling microscope (STM). 
In addition to STM the vortex core texture induced by 
the staggered chemical potential can be detected via a 
measurement of the orbital magnetic moment. This is 
because in this case the two core states carry opposite 
Q3 and it has been shown that in the presence of stag- 
gered potential the graphene quasiparticles carry a finite 
magnetic orbital moment which has opposite direction 
for the two valleys The complimentary nature of 
the texture of the two modes is important for differen- 
tiating the effect we are studying from any other such 
effect generated simply by translation symmetry break- 
ing. Also notice that changing the direction of magnetic 
field in the vortex core simply change the pattern to its 
complimentary form which is another unique feature of 
our result. 

In Fig. ||, we demonstrate the QVHE predicted by Eq. 
(||) . A single vortex is created at the center of the central 
SC graphene region and a finite supercurrent is passed 
along the strip by spatially twisting the phase of SC or- 
der parameter. Due to the interaction with the edges, 
the vortex texture is of the Kekule type. The asymmetric 
Kekule texture can be understood as the Q 5 accumula- 



FIG. 3: CDW texture of the vortex core states for (a) positive 
and (b) negative energy in the presence of a small staggered 
chemical potential around the core. Blue (red) points show 
surplus (deficit) quasiparticle density, with their size propor- 
tional to the magnitude. 




FIG. 4: The Kekule texture in the presence of supercurrent in 
y-direction. Blue (red) bonds show strengthened (weakened) 
bonds, respectively. The filled circle denotes the position of 
the vortex. 



tion induced by ramping up the supercurrent. This effect 
is a manifestation of the QVHE when one differentiates 
Eq. (0) with respect to Af . 

Although the accumulation of valley spin by supercur- 
rent can be viewed as a close relative of the QSHE, unlike 
the spin-filtered gapless mode at the QSH edge, there is 
no stable edge mode in the QVHE in general. In the QSH 
edge, time-reversal symmetry acting on spin protects the 
gapless nature of the edge modes, i.e. in order to make 
the edge modes gapped, one needs to break time-reversal 
symmetry, say, by adding magnetic impurities or induc- 
ing magnetic order at the edge. On the other hand, there 
is no such discrete symmetry acting on valley spin in the 
QVHE. In other word, quite generally, the edge mode of 
the QVHE is gapped because of valley spin symmetry 
breaking, leading to the Kekule or CDW texture at the 
edge. 

The connection between the physics discussed so far 
and the Landau forbidden transition discussed in Rcf. Q 
is the following. Imagine one starts with the supercon- 
ducting graphene and allows the vortices to proliferate. 
In the presence of proper symmetry breaking, which, say, 
favors the Kekule texture in the vortex core, such prolif- 
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oration will drive the SC into an insulator with Kckulc 
order. Since the symmetry groups of the two phases do 
not have the subgroup relation, the continuous transition 
triggered by the vortex proliferation defies the Landau 
rule. Conversely, if somehow the Kekule order is stabi- 
lized in graphene, as shown in Ref. 0, the vortex of the 
Kekule order can carry charge. As a result when the 
Kekule vortices proliferate and melt the Kekule order, a 
superconducting state emerges. The same scenario has 
been proposed originally for the putative direct transi- 
tion between VBS and Neel orders || [Hj]. The dual- 
ity between SC and VBS (Kekule) orders discussed in 
this paper can be mapped to that between the VBS and 
Necl orders upon a particle-hole transformation on one 
spin species. Other types of order parameters, which 
are seemingly different, and are not symmetry-related to 
each other in the Landau-Ginzberg sense, can be related 
in the same way; the quantum spin Hall phase (Kanc- 
Mele phase) and s-wave SC phase [^6j, Neel and <i-wave 
SC phases @, etc. 

The above duality is the physics of the Wess-Zumino- 
Wittcn (WZW) term. When several nearly degener- 
ate order parameters (or Higgs fields) compete to open 
mass gaps in a masslcss Dirac theory, one can derive a 
non-linear sigma model to describe the low energy/long 
wavelength order parameter fluctuations by integrating 
out the fermion fields. When the "mass term" asso- 
ciated with these orders obey certain anticommutation 
relation, a purely topological term, the WZW term, 
arises. For example, the five types of order parameter 
we have encountered (real and imaginary parts of su- 
perconductivity, two Kekule order parameters, and the 
CDW order) give rise to 5 mass-generating fermion bilin- 
ear j dr ^ , ^^4> a M a ^/ in graphene. The five hermitian 
mass matrices M„ anticommutc with each other, 



{M a ,M b } = 26 ab , a, 6=1,..., 5, 



(5) 



and also with the kinetic term, aoT x pop x + o z T y p z Py, of 
graphene. The WZW term generated upon integrating 
out the fermion fields fll§l is 



S 



WZW 



du / d*xe abcde 4> a d u 4> b d T( t> c d x $ d d y ^ e 



Area(5 4 ) Jo 

As discussed by Abanov and Wiegmann|l8| this term 



encodes the information on soliton charge, duality, etc. 
What is perhaps less appreciated is the fact that the 
physics of the WZW term is manifested even when the 
order parameters in question are not low energy degrees 
of freedom. Indeed, in the present paper, the SC or- 
der is artificially introduced by the proximity effect, and 
the Kekule/CDW order is not closed being realized. De- 
spite that the SC vortex core exhibits the texture of the 
dual order parameter. For graphene there are 36 pos- 
sible mass terms which anticommute with the graphene 
kinetic term, each representing a gapped order of some 
kind. Out of these 36 mass terms, one can form 56 sets 
of five mutually anticommuting mass terms Among 
these, SC-Kekule duality in this paper seems to be the 
easiest to probe experimentally. 
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